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Abstract 

We present a progress report for studies on maxima related to offspring in 
branching processes. We summarize and discuss the findings on the subject that 
appeared in the last ten years. Some of the results are refined and illustrated with 
new examples. 

1 Introduction 

There is a significant amount of research in the theory of branching processes devoted to 
extreme value problems concerning different population characteristics. The history of 
such studies goes back to the works in 50-ies by Zolotarev and Urbanik (see also 
0) who considered the maximum generation size. Our goal here is to summarize and 
discuss results on maxima related to the offspring. Papers directly addressing this area 
of study have begun to appear in the last ten years (though see "hero mothers" example 
in EJ.) 

Let A4 n denote the maximum offspring size of all individuals living in the (n — l)-st 
generation of a branching process. This is a maximum of random number of independent 
and identically distributed (i.i.d.) integer-valued random variables, where the random 
index is the population size of the process. M. n has two characteristic features: (i) the 
i.i.d. random variables are integer-valued and (ii) the distribution of the random index 
is connected to the distribution of the terms involved through the branching mechanism. 
These two characteristics distinguish the subject matter maxima among those studied in 
the general extreme value theory. 

The study of the sequence {A1 n } might be motivated in different ways. It provides 
a fertility measure characterizing the most prolific individual in one generation. It also 
measures the maximum litter (or family) size. In the branching tree context, it is the 
maximum degree of a vertex. The asymptotic behavior of A4 n gives us some information 
about the influence of the largest families on the size and survival of the entire population. 
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The paper is organized as follows. Next section deals with maxima in simple branch- 
ing processes with or without immigration. In Section 3 we derive results about maxima 
of a triangular array of zero-inflated geometric variables. Later we apply these to branch- 
ing processes with varying geometric environments. Section 4 begins with limit theorems 
for the max-domain of attraction of bivariate geometric variables. Then we discuss one 
application to branching processes with promiscuous matting. The final section considers 
a different construction in which a random score (a continuous random variable) is as- 
sociated with each individual in a simple branching process. We present briefly limiting 
results for the score's order statistics. In the end of the section, we give an extension to 
two-type processes. 

2 Maximum family size in simple branching processes 

Define a Bienayme-Galton- Watson (BGW) branching process and its n-th generation 
maximum family size by Zq — 1; 

Z n -1 

Z n = XAn) and M n = max XAn) (n = 1,2, ...), 

^-^ Kt<Z._i 
i=l 

respectively, where the offspring variables -Xj(n) are i.i.d. nonnegative and integer- valued. 

Along with the BGW process {Z n }, we consider the process with immigration {Z^ 71 } 
and its offspring maximum 

Z l n m = jr Xi(n) + Y n and M™ = mrc, M") (n = 1, 2, . . .), 

i—l ~ ~~ 

respectively, where {Y n , n — 1,2,...} are independent of the offspring variables, i.i.d. 
and integer-valued non-negative random variables. 

Finally, let us modify the immigration component such that immigrants may enter 
the n-th generation only if the (n — l)-st generation size is zero. Thus, we have the 
Foster-Pakes process and its offspring maximum 

Zl = V Xi(n)+I l7 o _ nX Y n and M° n = max X^n) (n=l,2,...), 

f l^n-l — u / Ki<Z° , 

i=l — — 

where I a stands for the indicator of A. 

Denote by F(x) — P(Xi(n) < x) the common distribution function of the offspring 
variables with mean < m < oo and variance < a 2 < oo. In this section, we deal with 
the subcritical (m < 1), critical (to = 1), and supercritical (to > 1) processes separately. 

2.1 Subcritical processes 

Let A4 n denote the maximum family size in all three processes defined above: {Z n }, 
{Z^ 71 }, and {Z°}. Let g(s) be the immigration p.g.f.. Also, let A n — {Z n _\ > 0} for 
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processes without immigration, and A n be the certain event - otherwise. The following 
result is true. 

Theorem 1 If < m < 1, then for x > 

lim P{M n < x\An) = 7(F(x)) (1) 

and 

lim E(M n \A n ) = V[l - l(F(k))} (2) 

71 — >00 * * 

k=0 

where 

(i) in case of {Z n }, 7 is the unique p.g.f. solution of j(f(s)) = m"f(s) + 1 — m and 
IS)/ 1 holds if in addition, EXi(n)\og{\ + Xi{n)) < 00. 

(ii) in case of process {Z^ 71 }, 0) holds provided £Jlog(l + Y n ) < 00 and 7 is the 
unique p.g.f solution ofj(s) = g(s)~f(f(s)). 0) is true if in addition, EY n < 00. 

(Hi) in case of process {Z®} we assume that i51og(l + Y n ) < 00. Then 7(s) = 
1 - E^ol 1 - 9{fn(s))] (0 < s < 1) and 7 (0) = {1 + £~ [1 - sC/UO))]}" 1 - Also, 0) 
holds if in addition, EY n < 00. 

Example 1 Consider {Z n } with geometric offspring p.g.f /(s) = p/(l — qs), where 
1/2 < p = 1 — q < 1. Then to = q/p < 1 and it is not difficult to see that j(s) = 
(1 — m)s/(l — ms). Hence 

(v - q)(l - q k+1 ) 
lim P (M n < k I Z n _i > 0) = f-^. 

It can also be seen ('201) that 

m m 
< lim E(M n \Z n > 0) < 



1 — pm n— >oo 1 — to 

Example 2 Consider {Z™} (see ioi#i 

/(s) = (1 + to-tos)" 1 (0<to<1) and g( s ) = f 1, (s) {v > 0). 
Then"f(s) = {{l — m)/{l — ms)Y, a negative binomial p.g.f ., and the above theorem yields 



1 - m 

1 — mF(j) 



< 



1 — to 



(1 — \ v 00 
— — ^ ) and lim KM™ = Vl- 
1 — mF{x) ) n—>oo ' ^—^ 

Example 3 Let /i = EY n . Consider {Z®} with 

f(s) = (1 + m — ms)^ 1 and g(s) = 1 — (/i/m) log(l + m — ms) (0 < m < 1). 
In this case 7 = (to, — /jlog(l — ms))/(m — /ilog(l — to)) and by the theorem 

m — /jlog(l — mF(x)) 



lim P{A1° < a;} 



TO — /i l0g(l — TO) 
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and 



to 



1 - m[(l + m) fe+1 



fc+11 



lim £W! = ,x *=° * < Mm 



m — fj, log(l — m) m — fj, log(l — to) 1 — to 

2.2 Critical processes 

In the rest of this section we need some asymptotic results for the maxima of i.i.d. 
random variables. Recall that a distribution function F(x) belongs to the max-domain 
of attraction of a distribution function H(x, 0) (i.e., F € D(H)) if and only if there exist 
sequences a(n) > and b(n) such that 

lim F n (a(n)x + b(n)) = H(x, 9) , (3) 

n — 'OO 

weakly. According to the classical Gnedenko's result, H(x;9) has the following (von 
Mises) form 

H(x;9) = exp{~h(x;9)} = exp{-(l + X0- 1 )- 8 } , 1 + x6~ l > 0; -oo < 6 < oo. (4) 

Necessary and sufficient conditions for F £ D(H) are well-known. In particular, F £ 
-D(exp{— x~ a }), a > if and only if for x > the following regularity condition on the 
tail probability holds 

1 - F(x) = x~ a L(x) , (5) 
where L{x) is a slowly varying at infinity function (s.v.f.). 

A. Processes without immigration. In case of a simple BGW process, the fol- 
lowing result holds. 



Theorem 2 Let to = 1 and a 2 < oo. (i) If 0) holds, then 

lim P [ Mn 7 \ {U) < x\Z n ^ >0)= _ (<i) 
'^oo \ a{n) ) 1 + a A h(x, 0)/2 

(ii) If 0) holds, then 

E(M \Z , ^ > 0) = _«ja > 
n-»oo ji 1 /"!! (n) sm(7r/o) ~ 

where L\{x) is certain s.v.f. with known asymptotics. 

The theorem implies that if F G £>(exp{— e~ x }) then the limiting distribution is 
logistic with c.d.f. (1 + e~ x ) ; and if F G Z?(exp{— x~ a }) then the limiting distribution 
is log-logistic with c.d.f. (1 + x~ a ) 1 . 
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Theorem 3 Let m = 1. a 1 = oo, and |3|) holds. Then for x > anc! 1 < a < 2 

1 



lim P ( . .. , M '\ < x | Z n _x > ) = 1 -j-, — -v, (8) 



which is a Burr Type XII distribution (e.g. '221) and 

V E(M n \Z n - 1 >0) _ 1 / 1 1 1 . 1 \ n . <9 , 

„i/[a(«-i)] i2 ( n ) a -l \a-l a(a-l)' a(a-l),/ 1 < ° ~ ; ' 

(9) 

where B(u,v) is the Beta function and L2(x) is certain s.v.f. with known asymptotics. 

Note that for a = 2 the right-hand sides in © (under assumption JUJ) and 
coincide with those in (0 and 0, respectively. The right-hand side in is the expected 
value of the limit in JBJ (see |2*2"|). 

Example 4 Let 1 — F (:r) ~ x _2 logx. /n £fe case one can check (see [201) that Theo- 
rem 3 with a — 2 implies 

lim P I 1/ ,,- M " >Q/0 < x I Z„-i > 0^ 4 ' 1 



71 — >OQ 



n^Qogn) 3 / 2 " ' / l + 4.x 



/or x > and 

E(M n \Z n -i >0) = 7T 
n 1/2 (logn) 3/2 2 



B. Processes with immigration {Z^}. Let /i = EY n . We have the following 
theorem. 

Theorem 4 Assume that m = 1, < cr 2 < oo, and < // < oo. (%) // 0) holds, then 

/M im -b(n) \ 1 
l im P -g^ < x = (10) 

1 o(n) " y (l + a 2 /i(z,0)/2) 2M/<7 



n — >oo 



J/ |0) is true, then 

lim - , " = ^ 4 + "!-- (o>2), 11 

where B(u,v) is the Beta function and Li2(x) is certain s.v.f. with known asymptotics. 

The theorem implies that if F € £>(exp{— e~ x }) then the limiting distribution is 

generalized logistic with c.d.f. (l + cr 2 e~ x /2) ; if F S £>(exp{— x~ a }) then the 

limiting distribution is a Burr Type III (e.g. [22]) with c.d.f. (l + a 2 x- a /2)~ 2/V<t2 . The 
right-hand side in (|llfl is the expected value of the limit in l|10|l (see |22p. 
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Theorem 5 Let m = 1, a 2 = oo, and holds. In addition, suppose 

Q(x):=- [ log[l - P(Zi m >0)}dt = c\ogx + d + e(x), (12) 

where lim.,.^00 e{x) — 0, c > 0, and d are constants. Then for x > 0, 

( M irn \ 1 
lim P[ , .. ; " <:r )=- r — - 1< a < 2), (13) 

which is a Burr Type III distribution (e.g. '221) and 

EJ\/V' 7rL f 1 1 \ 

lim — r7rl — t£ —=cB[c+— -.1 ; -] (1 < a < 2), (14) 

^^yji/IaM)]^^ ^ a(a-l) a(a-l)/ v ~ ; v 7 

where B(u,v) is the Beta function and L^{x) is certain s.v.f. with known asymptotics. 
The right-hand side in {1$ is the expected value of the limit in Wfy . 



Note that for c = 1 and a = 2 the right-hand sides in l|13(l and i|14|) coincide with those 
in (JSJ) and respectively. The condition l|12|) holds even when the immigration mean 
is not finite. Next example illustrates this point. 

Example 5 Following ]161j . we consider offspring and immigrants generated by 

f(s) = 1 - (1 - s)(l + (o - 1)(1 - 8 ))-V(«-i) and g(s) = exp{-A(l - s)^ 1 }, 
respectively. Then 0) holds and \12\) yields 

Q(t) = (A/(o - l))pog* + log(a - 1) + log(l + (a - l)*)" 1 ]. 

Therefore, 

( M lm \ 1 

1™ -P i /r ? Tti < X) = — -, tt-ytt — (1 < a < 2), 

n->oo \ n 1/[a(a_1)1 i3 (n) ~ J (1 + a;~ a ( a " 1 )) A /( a " 1 ) ~ 

and 

lim : /r f ; r = —B ( — J—,l-—±—) (l< a <2), 

where B(u,v) is the Beta function and L^{x) is certain s.v.f. with known asymptotics. 



C. Foster-Pakes processes {Z®}. The following limit theorem for under a 
non-linear normalization holds. 
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Theorem 6 Assume that m = 1, < a 2 < oo, and < /i < oo. If 

P(X 1 (l)>n) , . 

lim , -r = 1 (15) 

«^ooP(Xi(l) >n + l) v ; 

then for < x < 1, 

]imP (W(M°) <x \ (16) 
n— >oo y logn / 

where U(y) = 1/(1 -F(y)). 

Note that (|15|l is a necessary condition for X\ (n) to be in a max-domain of attraction. 

2.3 Supercritical processes 

Denote by M n (as in the subcritical case above) the maximum family size in all three 
processes: {Z n }, {Z*™}, and {Z^}. The following result is true. 

Theorem 7 Assume that m > 1 and EXi(n) log(l + JQ(n)) < oo. If |3J) holds, then 
n— »oo \ a(m") y 



Um -n/af^l 71 -n/a\ = I 1 ~ 

n->oo m I L\ [m n l a \ 



where L\(x) is certain s.v.f. with known asymptotics. 

(i) in case of {Z n }, ip is the unique, among the Laplace transforms, solution of 

^{u) = f^ium- 1 )), (u>0). (17) 

(ii) in case of {Z™ 1 } 7 we assume in addition that E\og(l + Y n ) < oo and 

oo 

V-M = II 9Mum- k )) (u > 0) , 
fc=i 

where (p(u) is the unique, among the Laplace transforms, solution of [17}) . 
(Hi) in case of {Z®}, we assume in addition that EY n < oo and 

oo 

1>(u) = g(<p(u)) - J> - }{^{um- n ))]P{Z^ =0) (u> 0) 
and (p(u) is the unique, among the Laplace transforms, solution of [17}) . 
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It is interesting to compare the limiting behavior of the maximum family size in the 
processes allowing immigration with that when the processes evolve in "isolation", i.e., 
without immigration. In the supercritical case, as might be expected, the immigration 
has little effect on the asymptotics of the maximum family size. The limits differ only 
in the form of the Laplace transform tp{u). In the subcritical and critical cases the 
mechanism of immigration eliminates the conditioning on non-extinction. Theorem 6 
for the Foster-Pakes process differs from the rest of the results by the non-linear norming 
of A4 n . The study of the limiting behavior of the expectation in this case needs additional 
efforts. 

It is known that some of the most popular discrete distributions, like geometric and 
Poisson, do not belong to any max-domain of attraction. This restricts the applicability of 
the results in the critical and supercritical cases above. A general construction of discrete 
distributions attracted in a max-domain is given in Wilms (1994). As it is proved there, if 
X is attracted by a Gumbel or Frechet distributions, then the same holds for the integer 
part [X] . Next we follow a different approach considering triangular arrays of geometric 
variables which leads to branching processes with varying environments. 

The results in this section are published in [T2]. and HH]-|231- In E3 an extension 
for order statistics is considered. 

3 Maximum family size in processes with varying en- 
vironments 

It is well-known that the geometric law is not attracted to any max-stable law. Therefore, 
the limit theorems for maxima in the critical and supercritical cases above do not apply 
to geometric offspring. In this section we utilize a triangular array of zero-modified 
geometric (ZMG) offspring distributions, instead. 

3.1 Maxima of arrays of zero-modified geometric variables 

In this subsection we prove limit theorems for maximum of ZMG with p.m.f. 



P{X i {n)=j) = ^ 




(n = 1,2,...) 



For a positive integer v n consider the triangular array of variables 



Ai(l),A 2 (l) 
Ai(2),A 2 (2) 



^(2) 



X 1 (n),X 2 (n), 



We prove limit theorems as v, 
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oo for the row maxima 



A4 n = max Xi(n). 



Ki<v, 
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Let A has the standard Gumbel law with c.d.f. exp(— e~ x ) for — oo < x < oo. 
Theorem 8 Assume that for some real c 

lim p n — and lim p n \og(v n a n ) = 2c. 

n — s-oo n — >oo 

A. If limji—yoo log(y n a n ) = oo, then c > and 

PnM n - \og{v n a„) -> A - c. 

B. 7/lirrin^oo log(^„a„) = a, (— oo < a < oo), then 

Pn M n ^(A + a)+. 

The idea of the proof is to exploit: (i) the exponential approximation to the zero- 
modified geometric law when its mean a n /p n is large; (ii) the fact that exponential law 
is attracted by Gumbel distribution. 

3.2 Processes with varying geometric environments 

Consider a branching process with ZMG offspring law defined over the triangular array 
above. Thus, we have a simple branching process with geometric varying environments. 
For this process we prove limit theorems for the offspring maxima in all three classes: 
sub-critical, critical, and supercritical. Define fiQ = 1, 

n 

p, n = E(Z n \Z = 1) = JJ nij (n > 1). 

i=i 

If the environments are weakly varying, i.e., p = linin^oo p n exists, then the processes 
can be classify (see |TT]) as follows. 

{supercritical if fj, = oo i.e. 5Z n (m n — 1) — > oo 

critical if fj, € (0, oo) i.e. ^2 n {m n — 1) < oo 

subcritical if ju = i.e. J2 n ( mn ~ -*-) — > ~ 00 

Define the maximum family size for the process with varying geometric environments as 
M 9 ° = max XAn), (n = 1, 2, . . .) 

l<i<Z„ 

In the result below the role played by v n before is played by B n _i where 



B n — ^ 



j=1 ^ 



Let V be a standard logistic random variable with c.d.f. (1 + e a: ) for — oo < x < oo 
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Theorem 9 Suppose that linin^oo B n — oo and for c real 

lim p n — and lim p n log(5 ra _ia„) = 2c. 

n — >oo n — >oo 

A. 7/lim„_ 00 log(B„_ 1 a„) = oo, then 

( Pn M 9 n e - log(B„_ 1 a„)|Z„_ 1 > 0) 4 V - c. 
5. // linin^oo log(-B n _ia n ) = a, (— oo < a < oo), then 

( Pn M n e \Z n ^ >0) 4 (V + a)+. 

Referring to the above theorem, we can say that the branching mechanism transforms 
Gumbel to logistic distribution. It is interesting to notice that this is in parallel with 
results for maximum of i.i.d. random variables with random geometrically distributed 
index discussed in [S]. 

Example 6 Let us sample a linear birth and death process (Bt) at irregular times. Let 
Z n = Bt n where < t n < i n +i — » <co < oo. If X and fi are the birth and death rates, 
respectively, and d n = t n — i n _i, then a n = m n p n , 



Pn = \ ' ^ m n =e" 



1 +\d n tf 1 -'*' 



and 



A( ^_~ 1} ifX^ 



{ Xt n if A = fl, 



I 1 ' /,._ = e ( A -M)t 



A. If X > fj, and 



lim — ^- = G [0, oo), 

ji^oo m n A — fi 



then 

( M 9e \ d 

[ iis_ - A - M t«) Vi>oUv-c. 

B. If X — fi and t n = n s l(n) (6 > 1), t/ien 

,A^R- lo ^' z "- 1>0 ) 4v - 

The results in this section can be found in 

4 Maxima in bisexual processes 

In this section we consider maxima of triangular arrays of bivariate geometric random 
vectors. The obtained results are applied to a class of bisexual branching processes. 
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4.1 Max-domain of attraction of bivariate geometric arrays 

The following construction is due to Marshall and Olkin (Hj. Consider a random vector 
(U, V) having Bernoulli marginals, i.e., it takes on four possible values (0,0), (0,1), (1,0), 
and (1,1) with probabilities poo, Poi, Pio, and pn, respectively. Thus the marginal 
probabilities for U and V are 

P(U = 0) =p + =Pao +Poi, P(U = 1) =Pi+ =Pw +Pn 
P(V = 0) =p+o =Pao +Pw, P(V = 1) =p+i =Poi + Px\- 

Consider a sequence {(£/„, V n )}^ =1 of independent and identically distributed with (U, V) 
random vectors. Let £ and r\ be the number of zeros preceding the first 1 in the sequences 
{Un}^ =1 and {Knj^i, respectively. Both £ and n follow a geometric distribution and, 
in general, they are dependent variables. The vector (£,77) has a bivariate geometric 
distribution with probability mass function for integer I and k 

[ PooPioP+o' _1 P+i if < I < k, 

p(t = i, v = k) = I p l ooPll _ if / = *, (is) 

{ PooPoiPo+ k ~ 1 Pi+ if < k < I. 

and 

P(£ > I n > k) - i P ™ lp +° l if ° " l ~ K (19) 

The marginals of £ and r\ for integer I and k are P(£ = I) = pi+p l 0+ (I > 0) and 
P(rj = k) = p + ip+ Q [k > 0), respectively and 

F e (0 = P(£>0=Po+ 1 C>0), F I) (A:)=F(r 7 >fc)=p^ 1 (fe > 0). (20) 

For n = 1,2,..., let be a positive integer and {(£i{n), r]i(n)) : i = 1,2,..., t^ n } 
be a triangular array of independent random vectors with the same bivariate geometric 
distribution l|18|) where pij are replaced by Pij{n) = 0, 1) for n = 1,2, . . . That is, 

(£i(l),m(l)),(6(l),»72(l)), .... (^(1)^(1)) 
(£ 1 (2),7 ?1 (2)),(6(2),r 72 (2)), (^(2)^(2)) 

Below we prove a limit theorem as i/ n — ■> 00 for the bivariate row maximum 



(A^^,A1^) = ( max &(n), max %(n)J 



Theorem 10 Lei lhrin^oo i/„ = 00. If there are constants < a, 6, c < 00, suc/i £/ia£ 
lira pu (n) log z/ n = 2c lim ^ r log v n — a and lim ^ 01 ^ - log ^„ — b, (21) 
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then for x. y > 

lim P (pni^Mi - log^„ < x, pn(n)Ml - logv n < y) 

n — >oc 

= exp {-e- x -°- c - e- y - b ~ c + e- n"«{*.v}-«-*-4 . 

Proof Set x n = (x + \ogv n ) /pn(n) and y n = (y + log u n )/pu(n). 

P{Mi<x n ,Ml<y„) = (F(x n ,y n )) Vn 

= (l - Ft(x n ) - F n (y n ) + P(£i(n) > x n , ^(n) > y n ))" n . 

Let x < y and thus, x n < y n . Taking logarithm, expanding in Taylor series, and using 
(TTOj) and we obtain 

\ogP(Mi<x n ,Ml<Vn) (22) 
= ^ n log(l -F e (x n ) -F v (y n ) +P(£ l (n) > x n ,r] t {n) > y n )) 
= -v n {[Fs(x n ) +F n (y n ) - P(&(n) > x n ,r)i(n) > y n )](l + o(l))} 

: : - (v n Po+(n) [x " ]+1 + v nP+0 (n) ly " ]+1 - ^Poo W)H [y " FK1 ) (1 + o(l)). 

Write [x n ] = x n — {x n }, where < {x n } < 1 is the fractional part of x n . It is easily 
seen that \im n ^ 00 (p Q+ (n))^ +1 = lim n ^ oo (p 0+ (n)) Xn+1 ^^ x ^ = ]hxL n ^ oo (p +(n)) Xn as 
n — > oo. Furthermore, taking into account (|21[) . we have 

log(p n pol(n)) = logv n + X+ ° gi/ " log(l - pi+(n)) 



= logi/ n - p ^ I pii(n) +pxo(n) + gCpnfa) +Pio(")) 2 + 0(pf + (n)) 

/ -. . h u Pio(n) , (Pn(n) +p 10 (n)) 2 2 
= + o 1 ^ logi/„ logu n + Pl n 

= -x(l + o(l)) (^gj + ip u („)) logi/ n (l + o(l)) + 0(p?!(n)) 
— > — x — a — c . 
Therefore 

lim v nPo+ (n)^ +1 = e- x ~ a - c . (23) 

n — >oo 

Similarly we arrive at 
lim zvP+oW [2M]+1 = crV - b - c and lim v n p m (n)^ +1 ^ e -x-a-b-c 

n — *oo n — >oo 

Finally, 

i »«-x„/ \ (y-iogi/ n ) - (z-iogi/ n ) 

l gpff« " ( n = _^ log 1 - pn (n) - poi (n)) 
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" " '" 'pii(n) +poi(n) + ~(pn(n) +Poi(«)) 2 + 0(p+i(n)) 



Pn(n) 



•'■ - - '/ - (// ■' ) ( + + («)(1 + + 0( P 2 +1 (n)) 



Thus, 



lim p +0 (7i) fe " 1 " [;r " 1 = e x " f . (25) 



The assertion of the theorem for x < y follows from i|22[l - i|25|) . The case y < x is treated 
similarly. This completes the proof. 
In particular, if a = b = then 

lim P ( Pll (n)Mi - log^„ < x, Pll (n)Ml - log^ <y) = exp {_ e -^».»}-"} . 

Note that in this case the limit is proportional to the upper bound for the possible 
asymptotic distribution of a multivariate maximum given in 0], Theorem 5.4.1. 

For the componentwise maxima, applying Theorem 10, one can obtain the following 
limiting results. If pi+(n) log v n — > 2c% < 00, then 

lim P(pu(n)Mi -logv n < x) = exp {-e~ x ~ Cl } . 

n — ►ao 

If pi+(n) log v n — > 2c2 < oo, then 

lim P(p lx {n)Ml-\ogv n <y) = exp{-e^- C2 } . 

n — >oc 

4.2 Bisexual processes with varying geometric environments 

Consider the array of bivariate random vectors {(£i{n), Vi{ n )) '■ i = 1,2, n = 
0, 1, . . .}, which are independent with respect to both indexes. Let L : 1Z + x 1Z + — > 1Z + 
be a mating function. A bisexual process with varying environments is defined (see |14|1 
by the recurrence: Z — N > 0, 

(Zf +1 ,Z^ 1 ) = ^te(n), ??l (n)) 

and 

Define the mean growth rate per mating unit 
r n j =j~ 1 E(Z n+1 \Z n =j) (J = 1,2,...) and /i„ = r iu /i = 1 (n = l,2, ...) 



Z n+1 =L(^ 1 ,Z^. 1 ) (n = 0,l,...). 



n-l 



i=0 
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Lemma ( 14 ) If 

Ef 1 -? 1 ) (26) 



then 

Z n 

lim — = W a.s., 

n^oo fj, n 

where W is a nonnegative random variable with E(W) < oo. 

//, in addition, there exist constants A > and c > 1 such that 

n+j-i 

j[ r tl >Ac n 3 = 1,2,...; n = 0,l,... (27) 

and there exists a random variable X with £'(Xlog(l + X)) < oo such that for any u 

P(X <u)<P < ^ (n = 0,l,...), (28) 

then P(W > 0) > 0. 

Further on we assume that (£i(n) , rji(n)) are i.i.d. copies of the bivariate geometric 
vector (£,77) introduced above and that the mating is promiscuous, i.e., 

L(£(n), 77(77)) = £(n) min{l, 77(71)}. (29) 

Theorem 11 Let {Z„} &e a bisexual branching process with varying geometric environ- 
ments and mating function \2ty) . If 

00 00 

Y[p+o(j)Po+(j) ^ and ^p+i(n) < 00 , (30) 

j=l n=0 

then 

lim — = W a.s., (31) 
where W is a nonnegative random variable with E(W) < 00 and P{W > 0) > 0. 

Proof To prove the theorem it is sufficient to verify the assumptions l|26|) - (|28[l in the 
above lemma. First, we prove that (H| holds. Indeed, for j > 1 

jr nj = E(Z^ +1 mm{l,Z^' +1 }) (32) 
= EE(ZF +1 min{l,Z™ +1 }\Z™ +1 ) 
= (l-P(Z 7 f +1 = 0))EZF +1 

- (1 -J ( n )) jP0+{n) 
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where we have used that both Z^ +l and Z^ +1 are negative binomial with parameters 
(j,p+i{n)) and (j,p 1+ (ri)), respectively. Thus, 



i , „^Po+(n) _ Po+(n) 



r n = lim r„j = lim (1 - p J +1 (n)) — - = — . (33) 

Now, (I32[l and l|33|) imply 1 — r n \/r n = p+i(n), which along with (|3l)[l leads to l|2tj|) . 
Let us prove idHJ • Indeed, for k > 1 

OO 

P(L(£(n), r?(n)) = k) = ]T P(£(n) min{l, ^(n)} = k\ v (n) = j)P(t](n) = j) 

3=1 

OO 

3=1 

OO 

= Pi+(«)Po+( n )X] p + 1 ( ri ) p +oW 

= P+o(n)Pi+(n)PQ + (n) . 

Therefore, P(£(£(n), ?j(n))/r n i > w) = p{,"^" 1 ' +1 (n) and hence, similarly to (|23[1 . taking 
into account l|3U|l . we obtain 

logPf r ^ >") ~ ur nl logp +(n) 

jj +0 (rt)7j 0+ (n) 
= -u 7 \ Pi+ « ! + o 1 

— > — u 

Thus, limn^oo P(L(£(n), r){n)/r n i) > u) — e~ u , which implies 

Finally, to prove (|27[) . observe that (|32[1 implies for any j and n 

n+j — l n+3 — 1 n-j-j — 1 

n = n p+o(*)po+(i) n 

»=j »=3 »=3 

oo n+j—1 

> Ac™ , 

where A = Ili^i-P+oWPo+M > (provided that the product in l|3U|l is finite) and 
c = rrrmj^j-p^ (i) > 1 (pn(i) — > under (|30|I % 1. (|27|l also holds if the product in (|3L)fl is 
infinite. Now, referring to the above lemma we complete the proof of the theorem. 
Define offspring maxima in the bisexual process {Z n } by 

(Ml,M™) = ( max d(n), max ^(n) 
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Theorem 12 Assume that p, n — > oo and there are constants < a, 6, c < oo, such that 
lim pn(n) log/i n = 2c lira ^ 10 ^ | log/i K = a and lim [ . log /i w = b. (34) 

ri-*oo n-»oo pn(n) n-^oo pn[n) 

Also assume that 

oo oo 

Y[ 'P+oU)po+{j) ^0 and ^p +1 (n) < oo . (35) 

i=l n=0 

T/ien 

/>oo 

lim P (pn(n)«M£ - log^„ < x,pu(n)M% - log/i n <y)= (G(x,y)) z dP(W < z), 
where 

G(x, y) = exp < — e — e " +e x >. 

Proof Set x n = (x+\ogfj, n )/pii(n) and y n = (y + logfi n )/pu(n). Under assumption 
()34p. Theorem 11 implies 

P {Ml < x n ,M™ < y n | Z n =k) = (F(x n ,y n )) k - (36) 

Under (|35|l . Theorem 12 implies 

lim P( — <#) = P(W<x). (37) 



Therefore, by and ljT71 . 



P[M n < -—- ,M n < -— = > P(Z n = k) (F(x n ,y n )) 



k=0 



jrp(— = —) (F(x n , Vn )r^ 

(G(x,y)) z dP(W <z). 

Next example, adopted from |10|. shows that the various conditions in Theorem 13 
can be satisfied. 

Example 7 Let a > 1 and f3 > 1. Sei 

p n (n) = n- Q and p 01 (n) = p 10 (n) = rr (a+/3) (n > 2). 
is not difficult to see that with this choice ofpij(n) — 0, 1), we have 
log /i„ ~ an log n as n — ► oo 
and 6o£/i (" with a = b = c = 0) and \35)) are satisfied. 

The exposition in this section follows JUli extending some of the results there. 
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5 Maximum score 



In this section we assume that every individual in a Galton- Watson family tree has a 
continuous random characteristic which maximum is of interest. 

5.1 Maximum scores in Galton- Watson processes 

Let us go back to the simple BGW process and attach random scores to each individual 
in the family tree. More specifically associate with the j-th individual in the n-th 
generation a continuous random variable Yj(n). Arnold and Villasehor (1996) published 
the first paper studying the maxima individual scores ("heights"). Pakes (1998) proves 
more general results concerning the laws of offspring score order statistics. Quoting |17| . 
" these results provide examples of the behavior of extreme order statistics of observations 
from samples of random size." Define by M(u\ n the fc-th largest score within the n-th 
generation and by M(u\ n the fc-th largest among the random variables {Yi(n) : 1 < 
* < Z v ,0 < v < n}, i.e., the fc-th largest score up to and including the n-th generation. 
Pakes (1998) studies the limiting behavior of "near maxima", i.e., (upper) extreme order 
statistics M( k y n and M^ n when n — > oo and k remains fixed. The two general cases 
that arise are whether the law of Z n (or the total progeny T n — J^™=o %y), conditional 
on survival, do not require or do require, normalization to converge to non-degenerate 
limits. 

If no normalization is required then no particular restriction need to be placed on 
the score distribution function S, but the limit laws are rather complex mixtures of the 
laws of extreme order statistics. The principal result states that 

oo k—1 / -\ 

Hm P(M [khn < x\A n ) = ]T ]T ( 3 \ (1 - S(x)) i S^ i (x)g j , 

j=l i=0 ^ ' 

where it is assumed that the conditional law Q n of Z n given A n (A n includes non- 
extinction) converges to a discrete and non-defective limit Q and gj denote the masses 
attributed to j by Q. 

If normalization is required then one must assume that the score distribution function 
S is attracted to an extremal law, and then the limit laws are mixtures of the classical 
limiting laws of extreme order statistics, let us assume that there are positive constants 
C n t oo such that for the conditional law Q n we have Q n {xC n ) => N(x), where N(x) is a 
non-defective but possibly degenerate distribution function. Assume also that the score 
distribution function S is in the domain of attraction of on extremal law given by (@J). 
The general result in ^2] is 

n^oo \ a[C n ) ) 7=t l - h 

Example 8 Consider an immortal (i.e., P{X = 0) = 0) supercritical process with shifted 
geometric offspring law given by its p.g.f. f(s) = s/(l + m — ms) (m > 1), then (see 
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Pakes (1998)) becomes 



hmP (M { ^-KC n )^ x{ X =i ( h(x, 0) 



rwoo y a (C n ) ~ 1 J \l + h{x,9) 

Thus the limit has a generalized logistic law when the score law is attracted to Gumbel 
law, h(x,8) = e~ x ; and a Pareto-type law results when S is attracted to the Frechet law. 

Phatarford (see |T2|) has raised the question (in the context of horse racing), "What 
is the probability that the founder of a family tree is better than all its descendants?" 

oo 

The answer turns out to be E (T 1 ^ 1 ) , where T = Z n is the total number of individuals 

n=0 

in the family tree. More generally, if t« is the index of the generation up to the n-th 
which contains the largest score, Pakes (1998) proves that 

P(r n = k) = E V (fc = 0,l,...,n), 

as well as limit theorems for r n as n — > oo. 
This subsection is based on [2] and |17j . 

5.2 Maximum scores in two-type processes 

Let each individual in a two-type branching process be equipped with a non-negative 
continuous random variable - individual score. We present limit theorems for the max- 
imum individual score. Consider two independent sets of independent random vectors 
with integer nonnegative components 

{X 1 (n)} = {(X 1 1 J .(n),X 2 1 i (n))} and {X 2 (n)} = X^n))} (j>l;n>0). 

A two-type branching process {Z(n)} = {(Zi(n), Z 2 (n))} is defined as follows: Z(0) ^ 
a.s. and for n = 1, 2, . 



Zi(n-l) Z 2 (n-1) 

z 1 ( n )= Yl x h( n )+ E 4w. 

i=i 3=1 

Zi(n-l) Z 2 (n-1) 

Z 2 (n)= X h( n )+ E *2». 

3=1 3=1 

Here X l k -{n) refers to the number of offspring of type k produced by the j-th individual 
of type i. With the j-th individual of type i living in the n-th generation we associate 
a non-negative continuous random variable (ij(n), (i — 1,2) "score", say. Assume that 
the offspring of type 1 and type 2 have scores, which are independent and identically 
distributed within each type. Define the maximum score within the n-th generation by 

Mi = maxjA^ 1 , M c n 2 }, where M% = max CijW (i = l,2). 

l<j<Zi(n-l) 
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Note that this is maximum of random number, independent but non-identically dis- 
tributed random variables. Let Fi(x) = P(d < x) (i = 1,2) be the c.d.f.'s of the scores 
of type 1 and type 2 individuals, respectively. 

Assumption 1 (tail-equivalence) We assume that F\ and F2 are tail equivalent, i.e., they 
have the same right endpoint x$ and for some A > 

hm = A. 

x^xa 1 — F2(X) 

Assumption 2 (max-stability) Suppose F\ is in a max-domain of attraction, i.e., Q holds. 

We consider the critical branching process Z(n) with mean matrix M, which is 
positively regular and nonsingular. Let M has maximum eigenvalue 1 and associated 
right and left eigenvectors u = [u\, U2) and v = (v\, V2), normalized such that u ■ v = 1 
and u • 1 = 1. 

Theorem 13 Let {Z(?i)} be the above critical two-type branching process. If the off- 
spring variance IB < 00 and both Assumptions 1 and 2 hold, then 

Hm P ( ^kllpM < x]Z{n) + 0) = — 1 — , (39) 

00 y a[v\Bn) J 1 + n(x,0) + (V2/V1) n(cx + a, 0) 

where if —00 < 9 < 00 is fixed, then c — A 1 /' 61 ' and d = 0; if 9 —> ±00, then c — 1 and 
d = \nA. 

Proof. Since F\(x) and F2(x) are tail- equivalent, we have (see P-67) 
lim (F 2 (a(n)x + b(n))) n -» H(cx + d,9), 

n — >oo 

where the constants c and d are as in (|39() . On the other hand, it is well-known (see [3], 
p.191) that for x > and y > 

lim P f < x, < y\Z(n) + 0) = G(x, y), 

where the limiting distribution has Laplace transform 

i>{\, 11) = 1 (A > 0, y. > 0). (40) 

1 + A + [i 

Set x n = a{v\Bn)x + b{v\Bn), s n = k/v\Bn, and t n — l/v2Bn. Referring to the 
definition of both and process {Z(n)} we obtain 

00 

P(Mi<x n \Z n ^0) = ]T P(Z(n) = (k,l)\Z(n)^0)P(m a x{M c n \Mi 2 } < x n ) 

(k,l)=o 

E„ f Zi in) k Z 2 (n) I ,„, . , \ r „ / , „ , 

P -V = —5", = TIT Z ™ + i F ^nt ^ 

^ q \v\Bn v\Bn v 2 Bn V2Bn ) 
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* — ' \v\Bn v 2 Bn J 



{k,l)=0 
/>oo />oo 



ff(a;, 9) s H{cx + d, 9)( V2/vi)t dG(s, t) 



JO 

oc poo 



^0 



exp <{ -sh{x, 6) - t^-h(cx + d,6)\ <IG{s. I ! 



i'2 



1 + 0) + — h(cx + d,( 

where in the last formula we used the Laplace transform of G(u,v) given in (|40|l . The 
proof is complete. 

The two examples below illustrate the kind of limit laws that can be encountered. 
Example 9 Let F\ and F 2 be Pareto c. d.f. 's given for xt > 9{ > and c > by 



Fi(xi) = 1 - 



(i=l,2). 



Note that the two distributions share the same value of the parameter c. It is not difficult 
to check that the limit is log-logistic given by 



lim P 



Mi 



n->oo yO^viBn) 1 /' 



m <x\Z(n) ^0 



1+ 1 



V2 
Vl 



Example 10 Let F± and F% be logistic and exponential c.d.f. 's given by 

1 



Fi(xi) = l-e X1 (0 < X\ < 00) and F 2 (x 2 ) — 



1 + e~ x * 



(—00 < x 2 < 00), 



respectively. It is known that both are in the max-domain of attraction of H{x) = 
cxp{— cxp{— x}} and share (see ^jj, p.91) the same normalizing constants a(n) = 1 
and b(n) = Inn. This fact, after inspecting the proof of the theorem, allows us to by- 
pass the tail- equivalence assumption and obtain a logistic limiting distribution, i.e, for 
— 00 < x < 00 



lim P {M c n - log(uiBn) < x | Z(n) ^ 0} 



1+ 1 



V2 
Vl 



The results in this subsection are modifications of those in ^3] . 
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